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We introduce a modified divergence law for the energy-momentum tensor in the theory of uni- 
modular relativity. Consequently, an additional equation for the scalar curvature follows from the 
divergence of the field equations. The equations of motion are derived and the weak-field, low- 
velocity limit is investigated. It is found that the gravitational acceleration acquires a term that is 
proportional to the gradient of the mass density. We show that this term can provide the additional 
acceleration observed on astrophysical scales without the need for dark matter. 
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I. INTRODUCTION 

Einstein's theory of general relativity (GR) [1[ has en- 
joyed overwhelming success since its inception, predicting 
phenomena that have been observed in our solar system 
and beyond, including the perihelion precession of Mer- 
cury, relativistic effects in the Hulsc- Taylor binary pul- 
sar B1913+16, and the existence of black holes. While 
GR reduces to Newtonian gravity in the weak-field, low- 
velocity limit, and hence agrees with observation for a 
wide range of phenomena, anomalous accelerations have 
been observed for decades. First discovered by Zwicky in 
the 1930's [1 [H, velocities on the galactic scale are much 
larger than those predicted by GR when the source of the 
gravitational field is taken to be the observed visible mat- 
ter (see also, e.g., [1] and [|[). Zwicky postulated, and 
it is now generally accepted, that a considerable amount 
of non-visible matter must be present in the extragalac- 
tic regime in order to provide the additional accelera- 
tion required to maintain these excessive velocities. This 
non- visible matter is commonly called dark matter and is 
believed to resolve acceleration discrepancies observed in 
systems ranging from dwarf spheroidal galaxies with visi- 
ble masses ~ 10 7 M Q to clusters of galaxies with observed 
masses ~ 10 14 M Q . Furthermore, dark matter is believed 
to play a key role in structure formation of the universe 
and primordial nucleosynthesis, and is believed to sig- 
nificantly affect the anisotropy of the cosmic microwave 
background. Excellent reviews of the dark matter prob- 
lem are given in Refs. [1,0]. 

Despite thirty years of laboratory experiments and as- 
tronomical observation, dark matter has never been ob- 
served directly ||; its existence is only inferred indi- 
rectly due to its purported gravitational effects on vis- 
ible matter. Modifications of gravitational theory have 
been proposed % El El El El 

that may eliminate 
the need for dark matter, and perhaps the most well- 
known is the modified Newtonian dynamics (MOND) 
theory [H El, E?| • MOND is characterized by an accel- 
eration scale and predicts departures from a Newtonian 



force law in the extragalactic regime where dynamical 
accelerations are small. Recently, a relativistic gener- 
alization of MOND, Tensor- Vector-Scalar (TeVeS), was 
proposed [l8[ that resolves some of the earlier problems 
of the MOND theory. However, TeVeS has not been ex- 
perimentally confirmed. 

In the following, we show that the introduction of a 
modified divergence law for the energy-momentum ten- 
sor in the unimodular theory of gravity results in an ad- 
ditional field equation. In the weak-field, low-velocity 
limit, the equations of motion acquire a term that is char- 
acterized by gradients in mass density. It is shown that 
these gradients can provide the additional acceleration 
observed on the galactic scale without the need for dark 
matter. The field equations reduce to GR when these 
gradients are negligible. Note that modified divergence 
laws are also discussed in [19, 0, 01 . 



II. FIELD EQUATIONS 

According to the theory of unimodular relativity [22| , 
the metric tensor is reducible under the general coordi- 
nate category into two nontrivial geometric objects: g, 
the determinant of the metric tensor, and 7^, the rel- 
ative tensor g^v/ {yf—g) 1 ^ 2 of determinant —1. The de- 
terminant determines entirely the measure structure of 
spacetime, while the relative tensor alone determines the 
null-cone or causal structure. Unimodular relativity as- 
sumes a background measure field \J—g — a(x a ) so that: 



(1) 



This condition is satisifed with the method of Lagrange 
undetermined multipliers in the action integral: 



S = J [{R + 2nL r 



f—g-a)} d 4 



(2) 



where L m is the Lagrangian density for the matter fields, 
A is the Lagrange undetermined multiplier, and n = 
(see also [H [24], [H, HI])- The curvature scalar R = 



j^Rftv is constructed from the connections: 
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(3) 



2 



where a comma denotes partial derivative. Variation of 
A gives y/— g = a(x a ) and variation of the quantities g^ v 
yields the field equations: 



Taking the trace of Eq. © gives: 



X = --(R-kT), 



(4) 



(5) 



and substituting this back into the field equations pro- 
duces: 
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kT» v - i/eiy" ) . !()) 



Taking the covariant divergence of these field equations 
gives: 



KT^ = -lg^ {R-kT) v 



(7) 



where the contracted Bianchi identity {G^Y V = 0) was 
used. In the common formulations of unimodular rel- 
ativity [H, [23|, [13, l2a ]. the covariant divergence law is 
assumed: 



However, this is an unjustified assumption in unimodu- 
lar relativity, which is highlighted by the fact that the 
measure field is non-dynamical and therefore unimodu- 
lar relativity is fundamentally not invariant under arbi- 
trary four-dimensional coordinate transformations. Con- 
sequently, we cannot assume that the action for matter 
in unimodular relativity leads to the vanishing of the co- 
variant divergence of the stress-energy tensor. Of course 
one may write the unimodular action for matter in a 
diffeomorphism-invariant form; however, one cannot as- 
sume that in so doing this will lead to the standard di- 
vergence law of GR. 

Therefore, one must permit generalizations of the stan- 
dard divergence law in unimodular relativity since T^ u v = 
is an unjustified assumption. One way to proceed 
is to re-examine the physical definition of divergence 
in light of the metric decomposition ©. The covari- 
ant divergence, in turn, relies on the parallel transport 
(defined by the connections) of energy-momentum since 
the divergence calculation involves the sum of the fluxes 
of energy-momentum through infinitesimal hypersurfaces 
separated by infinitesimal space-time intervals. Thus, we 
begin by examining the impact of the metric decomposi- 
tion © on the connection ([3]) . Substituting Eq. © into 
Eq. © we find: 



1 a/3 ~ 1 a/3 t L a/3' 



(11) 



where 



so that Eqs. © and © give: 

d\ 

dx^ 



= 0. 



(9) 



Hence, the cosmological constant emerges as a constant 
of integration in the field equations. This formulation has 
the attractive property that the contribution of vacuum 
fluctuations automatically cancels on the right hand side 
of Eq. © 

The primary difference between unimodular relativity 
and standard GR is revealed by evaluating the covari- 
ant divergence of the stress-energy tensor in each theory. 
In standard GR, the covariant divergence of the stress- 
energy tensor is fixed from the field equations: 



f 



- -Rg 



(10) 



Taking the covariant divergence of each side gives = 
0. In unimodular relativity, as in standard GR, the co- 
variant divergence also follows from the field equations 
(cf. Eq. ©). However, in unimodular relativity, the 
divergence is a function of the scalar curvature R. Since 
the unimodular field equations are traceless, the scalar 
curvature R is undetermined, and therefore one may not 
evaluate Eq. ([7]) without an additional condition. Typ- 
ically, this additional condition is taken over from stan- 
dard GR by simply forcing the covariant diver gen ce to 
vanish in unimodular relativity as well [HI, [H, HI, [25| . 
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ha\,/3 + 7/3A,a ~ Jaf3,\] (12) 

+ 8%a, a - y^ a p(T,x\ , (13) 



and are the normalized minors of 7 M „. We see that 
the connections separate into two terms, one that de- 
pends on derivatives of 7^,,, and another that depends 
on derivatives of a(x). Hence, parallel transport may 
be decomposed into two fundamental entities: parallel 
transport due to derivatives in the null-cone structure 
of space-time defined by T'^a , and parallel transport due 
to derivatives in the measure structure of space-time de- 
fined by T^. Substituting Eqs. © and (JTTJ) into T^ v 
we obtain: 
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where 



1 1 QI/ 



(14) 



(15) 



We see that the covariant divergence also separates nat- 
urally into two entities, one that looks like a divergence 
with respect to the null-cone structure alone, and an- 
other that depends on derivatives of the measure field. 
Note that T'" v vanishes because y/—g. x — and = 

TVa 1 G ,fi 
1 a/3 + — ■ 

As pointed out above, one may not adopt the addi- 
tional condition T^" = in unimodular relativity since 
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diffeomorphism invariance has been relaxed and we do 
not know the transformation properties of the unimod- 
ular matter action a priori. However, we may use the 
substructure of the connections and the covariant diver- 
gence identified above to postulate a generalized diver- 
gence law that is consistent with the relaxation of the 
invariance group. In so doing, we will find a likely can- 
didate for the matter Lagrangian in unimodular relativ- 
ity (cf. Eqs. ((26)) and ((27)1 ). We proceed by general- 
izing the law of parallel transport of energy-momentum 
so that the the null-cone and measure contributions may 
be weighted independently. For example, consider a law 
of parallel transport of energy-momentum based only on 
the derivatives in the null-cone structure defined by the 
connection T'£p. This leads to a divergence law based on 
the vanishing of Eq. (fT5)l alone: 



T^ v = 



In this case Eq. §8§ becomes: 



rpfll/ _ -V 

a 



Q 1 

2 4 y 



(16) 



(17) 



More generally, we may consider parallel transport that is 
weighted arbitrarily by each term in Eq. (jTT)> . This leads 
to the following generalization of the divergence law: 



= Z*L [AT'"' + BTg^"] + CT a "T'£ v , 
a 



(18) 



where A, B, and C are arbitrary constants. This reduces 
to the standard covariant divergence law when A = B = 
C = 0. This is equivalent to Eq. (JT7D for A = 3/2, 
B = —1/4, and C = 0. This divergence law implies that 
the action for matter is not a scalar invariant. When 
A = C = this divergence law leaves the usual covariant 
law unchanged for a traceless energy-momentum tensor 
(e.g., the electromagnetic energy-momentum tensor). 

The modified divergence law (|18p can be used with Eq. 
((7)) to derive an equation for the scalar curvature: 

[AT' 1 " + BTg^ v \ + CT av T'£ u = -i^-(^-T 

(19) 

When T^ v — this gives R = constant (with — unde- 
termined). Substituting Eq. ([5]) into Eq. ((T9)) we find: 



9X_ _ 2 ^u_ r AnTflv + BkTg im + 2CT a TZ, (20) 
ox v a 



nv 



which is equivalent to Eq. © in the absence of mat- 
ter. We see that matter described by T^" can introduce 
a space-time dependence into the cosmological constant 
term without the need for dark energy or quintessence 

In addition, Eq. (|18p predicts mass creation and anni- 
hilation; this results from gradients in the measure and 
null-cone fields in regions where the energy-momentum 
tensor is non-zero. If we consider a (cosmological) model 



where the only non-zero component of T^ v is T 00 = c 2 p, 
where p is the proper density of matter, then Eq. (TPS)) 
yields (assuming a diagonal metric): 



(V 



G -fpc 2 [A-Bg™] 



Cpc 2 



-7°°7oo,o 



2ii 1 



Bpc'g 



Cpc 2 



-7 Too.i, 



(21) 



where i = 1 ... 3 is not summed. The coefficients A, B, 
and C determine the net change in the energy-momentum 
density as the universe evolves. 

The equations of motion of a free dust particle follow 
from the matter conservation law ()18|) . Substituting the 
stress-energy tensor for dust T^ v — p^u^u" into Eq. 1(15)). 
where - "' ' 
that satisfies u^u 



-j— is the four-velocity (with r proper time) 



(7 



—c 2 we obtain: 
- [Au"u v 



Bc 2 g^] + Cu a u v T'^. 

(22) 

Contracting with and noting u^u^^ = we find: 



( A + B) u -^u v 



C 



u p u a u"Yt- (23) 



We see that the momentum density pou* 1 is conserved if 
A=-B and C = 0. Substituting Eq. (J23J) back into Eq. 
gives: 



l"u% = —B-^- [u' l u" + cVI 

1 



(24) 



Consequently, the equations of motion for dust are 

'dx^ dx" 



d 2 x< 1 . dx a dx 

a!i dr dr 



dT 2 



B- 



dr dr 



+ c 2 g^ 



C 



dx a dx v 
dr dr 



r'" + — v' 13 dxp dxfl 

a " c 2 a " dr dr 



(25) 



We see that an additional force, proportional to the 
derivatives of the measure and null-cone fields, modifies 
the equations of motion of a free dust particle. The equa- 
tions of motion depend on the coefficients B and C ex- 
plicitly and the coefficient A implicitly through Eq. p9[) . 
If B = C = then they reduce to the usual geodesic 
equations of motion. For the case C = the equations 
of motion follow from either the Lagrangian: 



1 



dx^ dx" 1 



(26) 



which is numerically equal to zero, or the Lagrangian: 

(27) 



L 



u B \ i g^v 



dx^ dx v 
dr dr 
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III. NEWTONIAN LIMIT 

We investigate the Newtonian limit by considering the 
motion of a free dust particle in a circular orbit in the 
field of a spherically symmetric, static, localized mass of 
total gravitational mass M. In addition, we assume that 
a spherically symmetric, static, matter density p(r) is dis- 
tributed in the space outside the localized source so that 
T 00 = c 2 p(r) (with T^ v = for all other components). 
We assume that the matter density outside the source M 
does not contribute appreciably to the traditional grav- 
itational force. Thus, M p <C M, where M p is the total 
integrated mass contribution from p(r) as observed from 
infinity. This scenario will serve as a model for the mo- 
tion of gas clouds in the disks of spiral galaxies at radii 
well beyond the visible disk's edge. 

In order to determine the radial force on a dust particle 
in a circular orbit we need to solve for under the 
assumptions stated above. Assuming a ~ 1, so that a — 
1 + e(r), where e(r) <C 1 we may write g^ v ~ <$Jf , where 

g$P is the solution of the field equations (|6|) due to the 
localized mass M. Since Eq. (J6j> in the absence of matter 
is equivalent to the free-field Einstein field equations with 
an arbitrary cosmological constant, gffi is: 



2m A , , ., 

1 -r 2 } c z dt 2 + 

r 3 



dr 2 



sin 



(28) 



■| is the cosmological con- 



where m = MG/c 2 and A = 
stant. 

Using Eqs. (p|l. ([25]). and 
radial acceleration of gravity in the limit 2p- <C c: 



28l) we can determine the 



d 2 x 1 
dr 2 



mc 



Ac 2 



a 



Cc 2 T 



on 



c 2 p,i 
4p 



(29) 



We see that the gravitational acceleration is modified by 
two terms, one proportional to the gradient of the mass 
density, and another proportional to the gradient of the 
scalar curvature. If we assume that the measure field is 
proportional to the matter density, then we may write 

(cf. Eq. mnn): 



Ri = qkTi = — anc 2 p^i, 



(30) 



or not the additional force is repulsive or attractive. For 
the case (a — 1) > this additional acceleration may 
be responsible for the acceleration discrepancies observed 
on astrophysical scales, thus removing the need for dark 
matter. Using Eq. (j3"Tj) , we can calculate the density 
distribution of visible matter that would produce a flat 
rotation curve beyond the visible edge of galaxies. The 
velocity for a bound circular orbit is: 



c 2 P,i 



= [ a -l) — ^-r 



(32) 



where we have ignored the term due to the cosmological 
constant. This will be constant as long as: 



{a- l)-r— r 
4 p 



mc 



1 o i 



(33) 



where i>o is the constant value of the rotation curve. Solv- 
ing the above differential equation yields: 



p(r) = p exp 



4m 



(a — l)r 



(34) 



where (3 
force 



4<-,1 



7T-T. In the limit that the gravitational 

is negligible, the density decays as r _/3 . As- 



suming (3 is on the order of unity and vq ~ 10 cm • s , 
we estimate (a — 1) ~ 10 -7 . 

In addition to the anomalous accelerations measured 
in the extragalactic regime, other evidence has been put 
forward in favor of dark matter. For example, image 
distortion due to gravitational lcnsing for a large number 
of distant background galaxies suggests that 90% of the 
matter of the foreground (lens) galaxies is invisible [28| . 
For the case C = we may use either Lagrangian PU)) 
or P7|) to obtain the equations of motion for massless 
particles: 



dp 2 



rA , dx a dx$ 
a @ dp dp 



B- 



dx^ dx v 
dp dp 



(35) 



where p is an affine parameter. Hence, in order to com- 
pute the Einstein deflection angle one needs to know the 
mass density profile (cf. Eq. (fll?|) ) through which the 
photons propagate. With the assumptions stated above, 
it is straightforward to show that the orbits of massless 
particles in a plane with polar coordinates {r = u -1 , cf)} 
around a point mass m satisfy (cf. [29[): 



where a is a constant. In other words, matter determines 
scale only to the extent that the gradient in the curvature 
scalar is proportional to the gradient in mass density. 
Consequently, Eq. ([2^]) becomes: 



mc 



Ac 2 lA c 2 p,i 

— r+{a- 1)- . 

3 4 p 



(31) 



Note that p x < for any reasonable matter distribu- 
tion, therefore, the sign of (a — I) determines whether 



d 2 u 



u — 3mu 



B da 
a du 



1 

k 2 



(36) 



PL 



where k = ^ 
case when 4 s - = 

du 



This is equivalent to the standard 



Therefore, mass determination via 
gravitational lensing will be inaccurate if the right hand 
side of Eq. (|3"6"|) is not included. 

As is well known, analyses of radio Doppler and rang- 
ing data from the Pioneer missions indicate that there is 
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an apparent anomalous acceleration do ~ 8 x 10 _8 cm/s 2 
directed towards the sun [3(3]. According to the conclu- 
sions above, this acceleration may be attributed to the 
gradient in the mass density. As before, we may cal- 
culate the density distribution that would produce the 
anomalous acceleration. Thus, we find: 

p(r)=poexp(-^^rJ, (37) 



where po is a constant. If (a — 1) ~ 10~ 7 as concluded 
above and r ~ 10 15 cm, then p would be roughly constant 
over this range of radii. On the other hand, assuming 
the term in the exponential is on the order of unity, then 
(a - 1) - 1CT 13 . 
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